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FOURIER  1  RANSFORM  AUTOCORRELATION 


INTRODUCTION 


Optical  Power  Spectrum  (OPS)  System.  ■  The  optical  Fourier  Transform  Auto¬ 
correlation  (FTA)  system  described  in  this  report  was  investigated  for  its  ability  to 
overcome  two  significant  problems  associated  with  the  usual  OPS  (Optical  Power 
Spectrum)  system.  These  problems  are  the  presence  of  phase  noise  in  the  object 
transparency,  and  the  diffraction  effects  arising  from  the  finite  sampling  aperture  within 
which  the  Fourier  transform  is  taken.  The  essential  difference  between  the  conventional 
OPS  system  and  the  FTA  system  is  that  the  OPS  system  detects  the  modulus  squared  of 
the  Fourier  transform  of  the  amplitude  transmitted  by  the  object  transparency,  and  the 
FTA  system  detects  the  Fourier  transform  of  the  modulus  squared  of  the  amplitude 
(the  intensity)  transmitted  by  the  object  transparency.  This  report  presents  an  analysis 
of  the  FTA  system  and  of  the  realizations  that  have  been  made  to  implement  the  theory. 
Experimental  data  is  lacking  because  the  large  amount  of  data  necessary  to  characterize 
a  transform  can  only  be  taken  with  an  automated  system. 

The  study  was  based  on  the  fact  that  the  Fourier  transform  of  the  intensity  power 
transmittance  of  an  object  transparency  is  equal  to  the  autocorrelation  of  the  Fourier 
transform  of  the  amplitude  transmittance  of  the  object  transparency.  This  can  be  ex¬ 
pressed  in  one  dimension  as1 


F'll^(x)P)  -  F[  (x)  ]  ®F[  (x)  ) 


(i) 


where  the  symbol  ®  denotes  correlation,  F  represents  a  Fourier  transform,  and 
\p  (x)  =  A  (x)  e'0(x)  (where  A(x)  and  <p(\)  are  real  and  i  =  y/-\)  is  the  complex  am¬ 
plitude  transmittance  of  the  object  transparency.  Also, 

t/t  (x)  i//*  (x)  =  li//(x)  I2  =  A2  (x)  (2) 


is  the  intensity  transmittance  of  the  object  transparency.  The  intensity  transmittance 
is  independent  of  any  phase  term  <A(x). 


1  R  .  B  racewcll,  Tht  Fourier  Transform  and  /is  Applications.  McGraw-Hill.  1965.  I  his  book  is  the 
reference  for  all  the  I  ouricr  transform  mathematics  used  in  (Ins  report. 
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The  Effect  of  Phase  Noise.  ■  In  this  section,  the  effect  that  phase  noise  and  a  finite 
sampling  aperture  have  on  the  detected  signal  in  an  OPS  system  is  described.  An  OPS 
system  diagram  is  shown  in  figure  1 . 
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Figure  1.  Schematic  of  a  Generalized  OPS  System. 


A  squaring  detector,  such  as  a  photodiode  array  or  photographic  film,  is  used  to  detect 
the  light  energy  falling  on  the  transform  plane.  A  squaring  detector  is  a  device  that 
senses  the  energy  flow  in  a  radiation  field.  The  name  comes  from  the  fact  that  the  energy 
density  is  given  by  the  square  of  the  field  amplitude.  This  means,  in  effect,  that  if  the 
Fourier  transform  of  the  object  function  g  (x)  is  given  by  F  [  g  (x)  ]  =  \jj  (0,  then 
the  detected  energy  is  given  by  \p  (0  ^  *  (0  =  A2  (f). 

If  the  object  function  contains  some  phase  noise  terms,  then  we  can  define  in  one 
dimension. 


g(x)  =  g,(x)  +  ig2(x)  (g,  (x)  and  g2(x)  real)  (3) 


where  g,  (x)  represents  the  amplitude  transmittance  of  the  transparency  resulting  from 
the  silver  grain  image  on  it  and  g2(x)  represents  the  phase  noise  caused  by  emulsion 
surface  relief  and  index  of  refraction  inhomogeneities  in  the  emulsion  and  film  base. 


Fourier  transforming  g  (x)  and  separating  the  real  and  imaginary  parts  of  the  transform 
yields 


Re  jF’lgfx)]}  =  l  gi  (x)  cos  (2trfx)  +  g2(x)  sin  (2irfx)  )  dx  (4) 

Im  \F  [g(x)||  =  i/^  [gx  (x)  sin(2rrfx)  +  g2(x)  cos  (2jrfx)  ]  dx  (5) 


Upon  detection,  only  the  real  terms  will  be  detected  and  these  will  be  squared.  By  de¬ 
fining  the  cosine  and  sine  transforms. 


c  [  g, (x)  1 

=  f°  g,  (x)  cos  (27rfx)  dx 

—  OO  1 

(6) 

S  1  g2(x)  ) 

=  g2(x)  sin  (27rfx)  dx 

(7) 

Then,  detection  will  give 

IRe  \F  (g(x)]  }  |2  =  C2[gl(x)  ]  +  S2[g2(x)  ]  +  2  C  [ g, (x)  ]  S[g2(x)](8) 

Expressing  equation  8  in  words;  the  odd  components  of  the  phase  noise  g2(x)  will  be 
transformed  into  an  amplitude  distribution,  which  will  be  added  to  the  even  components 
of  the  amplitude  due  to  gjfx)  and  will  be  mixed  with  the  desired  signal  from  g,(x) 
when  the  Fourier  transform  is  detected.  Equation  8  shows  that  phase  noise  is  a  serious 
problem  in  pattern  recognition  because  it  obscures  the  desired  information  due  to  gj  (x) 
in  the  detected  OPS. 
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The  Effect  of  the  Sampling  Aperture.  ■  The  aperture  used  to  define  the  area  of  a 
photograph  that  is  sampled  and  transformed  is  also  a  source  of  noise.  The  limited  region 
of  g(x)  means  that  the  Fourier  transform  is  no  longer  evaluated  over  all  space,  but 
rather  only  between  some  finite  limits  x,  and  x2.  This  limitation  causes  the  transform 
of  g(x)  to  be  convolved  with  the  transform  of  the  aperture  function.  This  aperture 
function  is  defined  somewhat  artifically,  so  that  we  may  still  have  a  Fourier  transform 
integral  from  to  Let  the  aperture  function  be  defined  as  centered  on  x  =  0,  so 
that 


x,  =  -x2  =  d. 


aperture 


j  1  for  |x|<dl 
1 0  for  |x|>d  ' 


(9) 


and  the  transform 


d 

u 


g  (x)  e'i27r^x  dx 


/  °°  rect  —  1  g  (x)  e'i2tffx  dx 
-oo  2d  J 


(10) 


Now,  we  can  make  use  of  the  fact  that 


F  [  rect 


2d  sine  (2df) 


(11) 


where 


sine  (x) 


sin  (7rx) 


(12) 


and  the  theorem  that  tor  any  two  functions  a(x)  and  b(x) 


F  [  a(x)  b(x))  =  F  [  a(x)l  */-'lb(x)l 


(13) 


where  the  asterisk  denotes  convolution. 
The  final  result  is 


oo 


rect 


g(x)e‘‘27r*x  dx 


2d  sine  (2df)  *  F  [  g  (x)  ] 


(14) 


In  an  OPS  system,  the  convolved  amplitudes  are  squared,  which  effectively  prevents 
the  removal  of  the  aperture  effect  from  the  detected  signal.  In  general,  deconvolution 
can  not  be  performed  without  retransforming  the  complete  complex  signal  into  a  multi¬ 
plication  problem,  after  which  the  aperture  function  can  be  divided  out.  However,  de- 
convolution  is  an  excessively  difficult  and  noisy  operation  and  is  not  done  in  practice. 
The  effect  of  the  aperture  function  is  most  evident  at  the  “d.c.  spike”  of  the  transform 
caused  by  the  average  or  zero  spatial  frequency  which  contains  most  of  the  energy  in 
the  transform.  The  transformed  aperture  function  around  the  d.c.  spike  is  removed  by 
subtracting  the  amplitude  field  owing  to  the  aperture.  However,  if  the  amplitude  field 
is  squared,  then  the  aperture  effect  is  multiplied  with  the  desired  signal  and  cannot  be 
removed  by  subtraction.  Because  the  aperture  function  cannot  be  deconvolved  or  even 
subtracted  in  the  zero  order,  much  effort  has  been  spent  on  apodization,  defined  as 
rolling  off  the  aperture  transmittance  as  a  function  of  radius  to  prevent  ringing  of  the 
aperture  transform.  Unfortunately,  apodization  also  removes  the  important  property 
of  translation  invariance  from  the  detected  transform  modulus  and  severely  restricts 
the  usefulness  of  the  OPS  system  as  a  pattern  recognition  device. 


The  Fourier  Transform  Autocorrelation  System.  ■  Figure  2  is  a  diagram  of  the 
FTA  system  designed  to  implement  equation  1.  The  system  is  essentially  a  Twyman- 
Green  interferometer,  whose  mirror  plane  is  imaged  by  a  telescope  into  the  object  plane. 
Light  transmitted  through  the  object  plane  is  then  collected  and  detected  by  a  photo¬ 
multiplier  tube  (PMT)  detector.  The  two  beams  of  the  interferometer,  one  from  each 
arm,  give  two  transforms  F  [  g(x)  ] ,  and  an  angle  between  the  beams  gives  the  correlation 
to  synthesize  the  function 


F  1  g(x)  1  *  F  l  g(x)  ] 


as  in  equation  1 . 


The  title  “Fourier  Transform  Autocorrelation”  is  somewhat  misleading  because, 
even  though  the  original  idea  for  this  system  came  from  equation  1 ,  it  has  been  found 
to  be  more  fruitful  to  analyze  the  FTA  system  as  a  detector  of  moire  patterns  formed 
between  a  projected  cosine  grating  generated  by  the  interferometer  and  similar  spatial 
frequencies  in  the  object  transparency.  Complicated  objects,  such  as  an  aerial  photo¬ 
graph,  contain  many  spatial  frequencies  and  directions  in  which  these  spatial  frequencies 
are  found.  The  determination  of  the  amplitudes,  phases,  and  directions  of  these  spatial 
frequencies  constitutes  taking  the  Fourier  transform  of  the  object  transparency.  The 
analysis  section  that  follows  will  show  how  the  FTA  device  takes  a  Fourier  transform 
and  will  discuss  details  of  implementation. 
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ANALYSIS 


The  analyses  contained  in  this  section  are  briefly  described  below  (All  of  the  des¬ 
criptions  refer  to  figure  2).  The  analysis  proceeds  from  the  most  general  and  mathe¬ 
matically  simplest  form  into  more  detailed  considerations  dictated  by  the  particular 
apparatus  used.  Analyses  1  and  2  are  concerned  with  the  heterodyne  methods  that  have 
been  tried.  Analyses  3  and  4  determine  the  system  spread  function,  (therefore  its  re¬ 
solution),  and  analysis  5  is  concerned  with  system  tolerances  and  mechanical  require¬ 
ments. 


Analysis  1 :  Wave  Amplitude  Analysis  is  a  description  of  the  case  where  mirror  M 1 
is  tilted  through  an  angle  6,  and  mirror  M2  is  linearly  translated  in  a  di¬ 
rection  normal  to  its  surface  at  a  constant  velocity  v. 


Analysis  2:  Effect  of  Piezoelectric  Mirror  Oscillation  shows  the  effect  of  oscillating 
mirror  M2  sinusoidally  rather  than  translating  it  with  constant  velocity. 
This  method  has  great  advantages  in  experimental  simplicity  and  in  ease 
of  analysis  of  the  detected  signal. 


Analysis  3:  Effect  of  Sampling  Aperture  on  Detector  Bandwidth  shows  the  effect  of 
the  system  aperture  on  the  bandwidth  of  spatial  frequency  detection  in 
the  Fourier  transform  plane.  This  in  effect  defines  the  system  point  spread 
function  in  the  spatial  frequency  or  radial  direction  in  the  transform  space. 


Analysis  4:  Object  Rotation  Effects  shows  the  effect  of  object  rotation  on  signal 
detection.  In  effect,  this  defines  the  system  point  spread  function  in  angle 
(a)  in  transform  space.  Together  with  analysis  3,  we  then  have  a  theoretical 
point  spread  function  for  the  system. 


Analysis  5:  System  Position  Resolution  Requirements  is  an  analysis  of  the  resolution 
needed  in  the  tilting  of  mirror  Ml  and  in  the  rotation  of  the  object  such 
that  a  Fourier  transform  resolution  consistent  with  the  device  spread  func¬ 
tion  may  be  obtained. 
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Wave  Amplitude  Analysis.  ■  The  FTA  device  initially  operated  by  tilting  mirror  M 1 
through  angle  9.  Mirror  M2  remained  normal  to  the  beam  but  was  linearly  translated 
along  the  optical  axis  at  a  velocity  v.  The  object  transparency  was  located  in  the  object 
plane  of  figure  2.  Angle  9  lies  in  the  plane  of  the  interferometer,  which  will  be  called 
the  x,  z  plane.  Thus,  z  is  defined  to  be  along  the  path  of  the  light,  x  is  positive  upward 
in  figure  2,  and  y  is  positive  out  of  the  paper  in  order  that  x,  y,  z  may  form  a  right- 
handed  coordinate  system.  Since  the  telescope  is  a  1:1  system,  an  angle  6  at  mirror  Ml 
will  be  imaged  into  the  same  angle  9  in  the  object  plane.  Mirror  M 1  is  imaged  into  the 
object  plane  by  the  telescope,  as  shown  in  dashed  lines,  so  that  a  tilt  9  does  not  result 
in  a  beam  translation  in  the  object  plane.  The  two  plane  wave  beams  from  Ml  and  M2 
give  the  following  normalized  amplitudes: 

Ml:  e'ik6x  where  k  =  and  X  =  wavelength  of  light  (632.8nm) 

A 

M2:  e'^vt  where  v  =  velocity  of  mirror  in  direction  of  light  path 
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Figure  2.  Fourier  Transform  Autocorrelation  System. 


(18) 


The  two  forms,  equation  1 7  and  18  show  different  properties  of  the  result. 

Equation  17  has  the  important  second  integral  in  the  form  of  a  correlation  of  g(x)  2 
with  cos  (kdx).  The  definition  of  correlation  is 


g®h  =  J^g(u)  h(u  +  x)  du 


(19) 


To  give  the  response  of  a  finite  system,  the  correlation  must  be  multiplied  by  the  system 
aperture  L.  This  action  does  not  invalidate  the  correlation,  but  only  causes  convolution 
of  the  result  by  a  spread  function  characteristic  of  th?  aperture  (see  analysis  3). 

Equation  18  shows  that  for  a  given  spatial  frequency  cos  (k0x)  there  are  two  de¬ 
tectable  time-varying  terms  which  appear  in  quadrature.  If  we  can  obtain  a  reference 
signal  of  the  same  temporal  frequency  as  cos  (kvt),  we  can  separate  the  last  two  integrals 
of  equation  1 6  by  phase  detection.  This  separation  is  very  important  since  the  form  of 
these  two  integrals  is  that  of  a  cosine  transform  and  a  sine  transform  of  lg(x)(2.  If  we 
can  separate  these  two  integrals  in  the  readout,  we  can  then  obtain  the  complete  Fourier 
transform  of  |g(x)|2  from  the  identity. 


F[h(x)]  =  C  [  E(x)  ]  +  i  S  [O  (x)  ]  (20) 


where  C  =  cosine  transform 
S  =  sine  transform 
E(x)  =  even  part  of  h(x) 


O(x)  =  odd  part  of  h(x) 


If  the  function  h(x)  is  completely  real,  as  |g(x)|2  is,  then  its  even  part  transforms  to  a 
real  function,  and  its  odd  part  transforms  to  an  imaginary  function.  Therefore,  the  even 
and  odd  parts  do  not  need  to  be  separated,  since  the  sine  and  cosine  transforms  will  do 
that  tor  us.  We  then  get  the  simplified  result 


F[  lg(x)|2  1  =  C[  |g(x)|2  1  +  i  S  [  lg(x)|2] 


(21) 


which  is  a  combination  of  the  last  two  integrals  of  equation  18.  To  examine  the  effect 
of  having  the  complete  Fourier  transform  of  |g(x)|2,  define  the  apertured  amplitude 
function  in  the  input  plane  as 


g  (x)  =  rect 


-^•j  [  A  +  gf(x)  +  igj(x)  ] 


(22) 


where  we  have  split  g(x)  into  its  average  amplitude  A,  its  real  varying  part  gr(x)  and 
its  imaginary  varying  part  gj(x).  n 

lg(x)|2  =  rect2  |^j[A  +  gr(x)  +  igi(x)][A  +  gr(x)_.gi(x)] 

2  (23) 

=  rect  |  $-)  [A2  +  g2  r(x)  +  g2  .  (x)  +  2Ag  r  (x)  ] 

Now  use  the  fact  that  rect2  j  X  |  =  rect  j  >  and  transform 

F  [  |g  (x)  |2  =  D  sinc(Df)  *  jA25(f)  +  F  [g2(x)J  +  F  [g?(x)]  +  2AF  [gp  (x)l  | 

=  Dsinc(Df)*  |A25(f)  +  F[|g(x)|2]  +  2AF[gr(x)]}  <24) 

The  FTA  apparatus  gives  this  Fourier  transform,  and  the  term  A2D  sine  (Df),  centered 
a  t  0,  can  be  digitally  subtracted.  This  process  removes  the  effect  of  the  aperture 
convolved  wdh  the  d.c.  spike  in  the  transform.  Although  this  is  not  equivalent  to  decon¬ 
volving  the  aperture  function  at  all  points  of  the  transform,  it  does  remove  the  most 
troublesome  aperture  effect  which  is  troublesome  because  of  the  extreme  intensity 
of  the  d.c.  term.  As  stated  earlier,  this  is  done  without  apodization,  and  the  position 
invariance  property  of  the  transform  modulus  is  retained. 
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Effect  of  Piezoelectric  Mirror  Oscillation.  ■  This  method  was  investigated  because  of 
experience  with  a  system  incorporating  the  idea  presented  in  analysis  1 ;  where  the 
mirror  M2  was  translated  at  a  constant  velocity  v  over  a  long  distance.  This  translation 
was  accomplished  by  moving  a  stage  with  a  1-inch  micrometer,  driven  by  a  synchronous 
60-RPM  motor.  Mirror  M2  was  mounted  atop  the  stage.  A  switch  mechanism  reversed 
the  motor  at  the  end  of  its  travel  in  each  direction,  and  data  was  taken  in  the  middle 
of  travel.  Two  problems  were  evident  with  this  system.  First,  the  stage  did  not  have 
straight  ways,  and  a  progressive  tilt  was  introduced  into  the  beam  as  the  stage  moved 
away  from  the  center  of  its  range  where  it  had  been  initially  adjusted.  Second,  the  syn¬ 
chronous  motor  did  not  give  a  constant  velocity  because  of  an  inherent  property  of  such 
motors  called  “cogging.”  The  only  way  to  eliminate  cogging  is  by  using  a  flywheel, 
which  would  be  difficult  to  fit  and  would  hamper  reversing.  These  problems  prompted 
the  following  analysis  in  which  mirror  M2  is  oscillated  by  a  piezoelectric  (PZ)  driver 
rather  than  linearly  translated  by  a  micrometer-driven  stage.2  A  diagram  of  the  elec¬ 
tronics  required  for  PZ  modulation  of  the  FTA  is  shown  in  figure  4. 

In  equation  15  replace  (v  t)  with 


A  sin  (2trft)  =  A  sin  (cot) 


where:  A  =  amplitude  of  oscillation  of  PZ 

f  =  frequency  of  oscillation  of  PZ 
to  =  27rf  =  phase  velocity 

Then,  by  a  similar  analysis,  we  get  the  detected  signal 

/lh(x)|2  dx  =  /|g(x)|2dx  +  coslkA  sin  (cot)l  /ig(x)|2  cos  (k0x)  dx 
L  L  L 

sin[kA  sin  (cot)]  /|g(x)|2  sin  (k0x)  dx 
I 


2This  idea  was  developed  from  lectures  in  Geometrical  Optics  by  Duncan  Moore  at  the  Institute  of  Optics  ot  the 
University  of  Rochester,  NY. 


The  amplitudes  of  the  cosine  and  sine  transforms  can  be  expressed  as  Bessel  series  ex 
pansions  using  the  following  identities. 


cos  IkA  sin  (cot))  =  J0(kA)  +  2J2(kA)  cos  (2cot)  +  2J4(kA)  cos  (4cot)  +.... 


sinlkA  sin  (cot)l  =  2J,  ( kA)  sin  (cut)  +  2J3(  kA)  sin  (3cot)  +  2J5(  kA)  sin  (Scot)  +.... 


This  is  a  very  useful  result  because  we  see  that  simple  frequency  bandpass  filtering  will 
give  the  coefficients  of  the  sine  and  cosine  transforms  without  a  phase  reference  signal 
or  phase  detection  systems.  More  precisely,  if  we  detect  the  fundamental  frequency 
(cot),  then  for  a  constant  amplitude  of  PZ  oscillation, 


detected  amplitude 
2J,  (kA) 


(26) 


Detection  of  the  first  harmonic  (2cot)  frequency  then  gives 


detected  amplitude 

— -  =  N 

2J2  (kA) 


(27) 


where  M  and  N  are  the  amplitudes  of  the  sine  and  cosine  transforms  respectively,  of 
the  particular  spatial  frequency  of  |g(x)|2  which  is  being  tested  for  as  determined  by 
the  mirror  angle  0.  The  value  of  J|  (kA)  and  J2  (kA)  must  be  accurately  known  to  get 
a  good  result.  It  is  useful  to  plot  J,  (kA)  and  J2  (kA)  to  find  the  region  of  kA  that 
will  give  the  best  result.  For  kA>3,  the  Jj  and  J2  functions  are  approximately  in  quad¬ 
rature;  so  it  is  useful  to  consider  the  region  l<kA^3,  which  encompasses  amplitudes  of 
slightly  less  then  14X  to  around  ViX  at  which  kA  =  * .  The  graph  shows  that  the  region 
about  haltway  between  !4X  and  V2X  gives  the  best  combination  of  high  amplitude  and 
low  slope  for  the  functions  (figure  3).  Low  slope  is  desirable  so  that  the  value  of  the 
function  is  not  too  sensitive  to  slight  amplitude  changes  of  the  PZ  driver. 

Without  using  electronic  phase  detection  it  is  difficult  to  get  an  accurate  value  for 
the  optical  path  difference  (OPD  =  kA).  Without  an  accurate  OPD  it  is  impossible  to 
calculate  accurate  values  of  J,  (kA)  and  J2  (kA)  and.  thereby  to  get  good  values  for 
the  amplitudes  of  the  sine  and  cosine  transforms. 


The  principal  problem  associated  with  OPD  detection  is  that  the  origin  of  coordinates 
for  the  sinusoidal  fringe  pattern  projected  by  the  interferometer  is  unknown.  This  leads 
to  a  phase  ambiguity  that  is  not  constant  for  different  spatial  frequencies.  If  the  modulus 
of  F  [  |g(x)|2  1  is  required,  this  phase  ambiguity  does  not  have  to  be  recovered:  but  if 
the  complete  transform  is  needed  (as  it  is  for  such  operations  as  deconvolution  of 
aperture  effects),  then  a  phase  reference  must  be  provided.  A  reference  phase  can  be 
derived  at  low  spatial  frequencies  by  sampling  a  portion  of  the  field  input  to  the  object 
transparency.  This  could  be  done  with  a  beamsplitter,  a  microscope  objective  to  magnify 
the  fringes,  and  a  detector  to  sample  the  field  in  a  small  region.  Although  high  spatial 
frequencies  cannot  be  phase-referenced  by  the  method  described,  it  is  only  necessary 
to  know  the  phase  of  low  spatial  frequencies  (up  to  about  10  cycles  per  millimeter) 
for  the  purpose  of  deconvolution  of  an  aperture  function.  The  problem  of  aperture 
function  deconvolution  will  be  addressed  in  the  next  section  of  this  report. 

1  7 


TO  DIGITIZER 
AND  COMPUTER 


Figure  4.  FTA  System  Electronics  Block  Diagram. 
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Effect  of  Sampling  Aperture  on  Detector  Bandwidth.  »  The  detected  moire  pattern 
will  cause  a  signal  amplitude  in  the  detector  which  is  proportional  to  the  original  moire* 
fringe  contrast,  fringe  spacing,  and  system  aperture.  The  system  aperture  is  the  sampling 
aperture  on  the  object  and  will  be  circular.  Let  the  moire  pattern  equation  be 


I  =  I0  [  1  +  m  cos  (2nfx)  ] 

m  =  contrast  (28) 

f  =  spatial  frequency  in  x-direction 


Normalize  the  intensity  by  I  =  1  and  integrate  this  function  over  a  circular  aperture 
of  radius  C/2  using  polar  coordinates  r,  6  where  x  =  r  cos  0.3 

1  =  2  f*  2  [  1  +  mcos(2jrfx)]  rdrdfl  (29) 

2  i  {  ^  r  dr  dd  +  2  m  £  f^2  cos  (2rrfr  cos  6  +  27rfx0)  r  dr  d0 

The  arbitrary  phase  shift  2n  fxQ  has  been  added  to  make  the  result  general.  The  first 
integral  is  just  twice  the  detector  area.  Split  the  second  integral  into  two  integrals  by 


cos  (A  +  B)  =  cos  A  cos  B  -  sin  A  sin  B 


The  resulting  integrals  must  be  evaluated  from  tables,  but  the  final  result  for  the  detected 
signal  is 


output  =  (constant)  [  1  +  m 


2J,  (jrGf) 


cos(27rfx0) 


rather  than  just  m,  which  shows 


The  observed  fringe  contrast  is  now  m  -=J-i  (---°  rather  than  just  m,  which  shows 

\  7rCf 

that  as  the  aperture  is  made  larger,  the  Af  range  over  which  a  good  response  is  seen  is 
reduced.  The  effect  is  that  the  instrument  spread  function  in  the  synthesized  Fourier 
transform  space  is  a  Bcsinc  function,  just  as  it  is  for  conventional  optical  Fourier  trans¬ 
forms  using  circular  apertures. 


3This  section  developed  from  lectures  in  Interferometry  by  M.  Parker  Givens  at  the  Institute  of  Optics  of  the  University 
of  Rochester,  NY. 


Since  we  are  dealing  with  a  Besinc  function,  Rayleigh’s  criterion  may  be  used  for 
resolvable  points  in  the  transform  plane.  This  criterion  is  defined  as  the  distance  in  fre¬ 
quency  space  from  the  central  maximum  of  the  Besinc  function  to  its  first  zero.  The 
first  zero  of  the  function  is  at  l.220ir.  Substituting  the  argument  of  our  derived  function, 
we  get 


TrCf  =  1.220tt 


Af 


The  following  table  was  generated  using  Rayleigh’s  criterion: 


aperture  diameter  ( C)  in  mm 
1 

Am 

3 

4 

5 

6 

7 

8 

9 

10 


t  Af  in  cycles/mm 

±  1.22 
±  0.61 
±0.41 
±  0.31 
±  0.24 
±  0.20 
±  0.17 
±  0.15 
±  0.14 
±  0.12 


This  table  indicates  the  extreme  sensitivity  of  the  system.  The  narrow  spread  function 
means  that  the  tilt  of  mirror  Ml  to  scan  in  frequency  must  be  very  precise,  and  the  ro¬ 
tation  of  the  object  transparency  to  scan  in  direction  a  must  be  similarly  precise.  An 
exact  calculation  of  these  requirements  is  in  Analysis  5. 
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Object  Rotation  Analysis.  ■  For  the  purpose  of  this  analysis,  assume  that  the  sam¬ 
pling  aperture  of  the  FTA  is  infinite.  This  assumption  is  equivalent  to  saying  that  the 
device  point-spread  function  (psf)  will  be  a  delta  function.  To  find  the  psf  of  a  real  device 
with  a  finite  sampling  aperture,  we  only  have  to  convolve  the  psf  found  in  this  analysis 
with  the  Fourier  transform  of  the  sampling  aperture.  In  the  previous  section  this  function 
was  shown  to  be  a  Besinc  for  circular  apertures. 

The  assumption  of  an  infinite  sampling  aperture  simplifies  this  analysis  and  enables 
the  result  to  be  more  general.  The  detector  acts  as  a  low-pass  spatial  frequency  device 
that  will  only  pass  a  zero  spatial  frequency  moire'  for  an  infinite  sampling  aperture.  This 
low-pass  filtering  means  that  the  difference  frequency  between  the  cosine  grating  pro¬ 
jected  by  the  interferometer  and  some  spatial  frequency  in  the  object  must  be  zero  to 
be  detected.  To  derive  the  conditions  under  which  this  difference  spatial  frequency  is 
zero,  assume  an  object  function  consisting  of  a  single  spatial  frequency  in  the  x  direction. 
The  field  immediately  behind  the  object  transparency  is  then  given  in  normalized  form  by 

i //( x)  =  [  1  +  A  cos  ( 2 rr  f,  x)  ]  cos  ( 2  7r  f 2 x  +7)  (31) 

where 


A  -  transparency  modulation  as  a  fraction  of  average  transmittance 
f(  =  object  spatial  frequency 
f2  =  spatial  frequency  projected  onto  object 
7  =  arbitrary  phase  offset 


For  simplicity  in  notation,  define 


P 

0 


=  27rf(x 
=  27rf2x 


The  optical  field  is  then  rewritten  as 

\p  (x)  =  [  1  +  Acos(P)]  cos  (Q  +  7) 

(32) 

=  cos  (Q  +  7)  +  A  cos  (P)  cos  (Q  +  7) 

Detection  of  this  field  then  gives  (33) 

li// (x)|2  =  cos2  (Q  +  7)  +  2A  cos  (P)  cos2  (Q  +  7)  +  A2  cos2  (P)  cos2  (0  +  7) 


Since  it  must  collect  all  the  light  passing  through  the  FTA  device,  the  detector  must 
also  be  infinite  in  extent.  This  requirement  is  of  course  relaxed  for  a  real  device  with 
a  finite  sampling  aperture.  Use  trigonometric  identities  to  remove  squared  cosines,  and 
get  the  following  equation  whose  terms  are  grouped  into  significant  categories: 


li/Tx)!2  =  1/2  +  A2/4 


bias  terms 


+  1/2  cos  (2Q  +  27)  +  Acos(P) 

+  (A2/4)  cos  (2Q  +  27)  +  (A/4)  cos  (2P) 


(34) 

high-frequency  terms 


+  A  cos  (P)  cos  (2Q  +  27) 

cross  terms 

+  (A2 /4)  cos  (2P)  cos  (20  +  27) 


The  detector  acts  as  a  low-pass  filter  for  spatial  frequencies  so  that  the  high  frequency 
terms  will  not  be  detected.  To  determine  detectability  of  the  cross  terms,  expand  them 
by  trigonometric  identities. 

(35) 

A  cos  (P)  cos  (2Q  +  27)  =  (A/2)  cos  (P  -  20  --27)  +  (A/2)  cos  (P  +  2Q  +  27) 

2/4)  cos  ( 2 P)  cos  (20  +  27)  =  (A2/8)cos(2P- 2Q- 27)  +  (A2/8)cos(2P  +  2Q  +  27) 

(36) 


Low-pass  filtering  will  not  detect  the  sum  terms.  The  detectable  terms  are  now 
reduced  to  a  manageable  size.  Substituting  back  for  P  and  0,  we  get  the  detectable  terms 


1  *  (x)|2fi.,ered 


2  +  A2 

— - -  +  (A/2)  cos  [  27rx  ( f,  —  2f2 )  —  2-y J 


(37) 


+  (A2/8)  cos  [47tx  (fj  —  f2)  -  2-y ] 


Assuming  an  infinite  sampling  aperture,  the  cosine  terms  will  be  detected  when 


f,  -  2f2  =  0 


and 


ft  -  f2  =  0 


These  conditions  define  two  delta  functions  in  the  Fourier  transform  and  give  the  unfor¬ 
tunate  result  that  the  FT  A  point-spread  function  is  double  valued. 


A  physical  interpretation  of  the  above  result  can  be  seen  if  the  two  projected  light 
beams  are  diffracted  into  plus  and  minus  orders  by  an  object  function  consisting  of  a 
biased  cosine  grating,  which  gives  the  two  coincidence  conditions  as  shown  in  Figure  5. 


a)  condition  f,  -  2f2  =  0  where  f, 


object  spatial  frequency 


t,  -  projected  grating  spatial 
frequency 


INTERFERING  ORDERS 
-1  OF  BEAM  1 
+1  OF  BEAM  2 


INTERFERING  ORDERS 
0  OF  BEAM  1 
+1  OF  BEAM  2 

INTERFERING  ORDERS 
-1  OF  BEAM  1 
0  OF  BEAM  2 


Note  that  beam  1  is  fixed  and  beam  2  is  steered. 

Figure  5.  Physical  Interpretation  of  Double-Valued  Output  of  FTA. 
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Notice  that  the  condition  f,  -  2f2  =0  appears  as  interference  along  a  single  direc¬ 
tion  and  the  condition  f,  -  f2  =0  appears  as  interference  along  two  directions.  The 
influence  of  the  fundamental  term  (f,  -  f2  =  0)  drops  as  the  modulation  of  the  object 
grating  decreases  as  can  be  seen  from  the  table  below. 


A_ 

A/2 

A2/8 

A2/8  as  v 

1 

0.5 

0.125 

25 

0.5 

0.25 

0.031 

13 

0.3 

0.15 

0.011 

7 

0.1 

0.05 

0.0013 

3 

0.05 

0.025 

0.00031 

1.3 

The  variable  A  is  the  object  modulation,  A/2  is  the  coefficient  of  the  first  harmonic  term 
(fj  -  2f2  =  0)  and  A2 /8  is  the  coefficient  of  the  fundamental  harmonic  term. 

An  attempt  to  change  the  FTA  from  a  double-value  to  a  single-valued  device  would 
necessitate  Fourier  transforming  the  output  and  aperturing  in  the  detector  space  to  block 
out  unwanted  orders.  Examination  of  figure  5  shows  that  the  -1  order  of  beam  1  takes 
part  in  both  condition  (a)  and  one  of  the  beams  of  condition  (b).  Detection  of  an  un¬ 
ambiguous  frequency  requires  an  attempt  to  isolate  the  beam  corresponding  to  the  zero 
order  of  beam  1  and  the  +1  order  of  beam  2  as  shown  in  condition  (b).  Isolation  of  this 
desired  beam  is  made  easier  by  the  fact  that  beam  1  is  fixed  in  direction  so  the  aperture 
need  not  be  moved  as  the  device  scans  in  spatial  frequency.  Other  considerations  mitigate 
against  such  a  solution  however.  First,  the  time-varying  amplitude  of  the  isolated  beam 
will  have  a  coefficient  of  A2/16  and  will  give  signal-to-noise  problems  at  the  low  modula¬ 
tions  found  in  aerial  photography.  Second,  the  zero  order  of  beam  1  is  very  strong 
(bright)  and  will  cause  a  large  bias  term  leading  to  an  increase  in  detector  noise.  These 
two  considerations  are  very  important  because  bench  work  has  shown  that  the  FTA 
device  has  a  signal-to-noise  problem  that  is  quite  severe. 

If  it  is  desirable  to  attempt  to  aperture  the  output  of  the  FTA  device  as  described 
above,  a  criterion  for  minimum  detectable  spatial  frequency  can  be  derived.  This  criterion 
will  necessarily  be  dependent  on  the  size  of  the  sampling  aperture  used  in  the  object 
plane  because  readings  of  spatial  frequency  will  be  degraded  by  energy  from  unwanted 


h  —  ■  •  ■  m  ■  ( 


diffraction  orders  entering  the  detector  aperture.  This  energy  comes  from  ringing  of  the 
transform  of  the  input  aperture  function.  As  an  example  of  such  calculation,  note  that 
the  Besinc  function  has  the  normalized  form  given  in  the  preceding  section. 


2 J ,  (TTfif) 

"w 


where  C  -  diameter  of  object  plane  aperture 
f  =  spatial  frequency 


If  we  want  to  collect  energy  from  this  “spread  function”  out  to  the  third  zero-crossing 
without  collecting  energy  beyond  the  third  zero-crossing  of  the  adjacent  (and  unwanted) 
order,  we  have  the  situation  shown  in  figure  6. 


Zeros  of  Besinc  functions 
are  indicated  by  circles,  the 
sampling  aperture  is  in¬ 
dicated  by  a  dashed  line. 


Figure  6.  Placing  an  Aperture  on  the  FT  A  For  Single-Valued  Output. 


From  tables,  the  third  zero  of  the  Besinc  function  is  found  when  the  argument  is  10. 1 73, 
that  is 


2  6 


7T  C  f 


10.173 


The  orders  are  separated  by  twice  this  distance;  so  the  solution  for  minimum  detectable 
spatial  frequency  is 

20.35 
f  =  - 


which  can  be  seen  to  be  inversely  related  to  the  object  sampling  aperture  size.  Some 
representative  values  are 


Notice  that  the  optimum  size  of  the  detector  aperture  is  also  dependent  on  the  object 
sampling  aperture  diameter.  Criteria  for  any  separation  of  desired  and  unwanted  orders 
can  be  similarly  derived. 


System  Position  Resolution  Requirements.  •  To  realize  the  potential  device  resolution 
given  by  Rayleigh’s  criterion,  mechanical  rigidity  and  high  positional  accuracy  of  the  tilt 
(0)  of  mirror  Ml  (figure  2)  and  the  rotation  (0)  of  the  object  transparency  in  its  own 
plane  will  be  required.  The  analysis  in  this  section  will  provide  some  numbers  to  be  used 
in  apparatus  design.  Assume  that  the  maximum  sampling  aperture  to  be  used  in  the 
system  will  be  10  millimeters  in  diameter,  that  the  maximum  spatial  frequency  of  interest 
is  100  cycles  per  millimeter,  and  that  the  telescope  in  the  device  images  at  1:1. 


A  maximum  aperture  of  10  millimeters  implies  by  Rayleigh’s  criterion  that  the 
frequency  spread  Af  =  ±0.122  cycle  per  millimeter.  This  then  corresponds  to  the 
necessity  of  setting  the  angle  of  mirror  M2  to  a  A 0  given  by  differentiating  the  grating 
equation  for  normal  incidence. 


sin  0  =  f  X  where  f 


1 _ 

grating  spacing 


d  (sin  0) 
df 


X 


cos  0  d 0 
df 


X 


(38) 


cos  0 

0.6328  x  10‘3 

d0  =  -  mm  (39) 

0.122 

d0  =  7.7202  x  10‘5  radian 

“  16  arc  seconds 


Using  a  micrometer  that  can  be  read  to  0.001  inch,  we  need  a  pivot  arm  on  mirror  Ml 
of  a  length  given  by 


d0  (radians) 


.001 

7.7202  x  1 0'5 


.001 

x 


x 


1 2.95  inches 


This  will  give  a  A0  least  count  consistent  with  Rayleigh’s  criterion.  Using  a  micrometer 
with  a  total  travel  of  1  inch  will  give  a  maximum  angle  of  about  4°,  and  since  3.63* 
corresponds  to  a  grating  frequency  of  100  cycles  per  millimeter  this  will  be  a  good  choice. 

A  standard  micrometer  screw  has  40  turns  per  inch.  If  the  micrometer  is  turned  by 
a  60  revolutions  per  minute  motor,  then  a  sample  must  be  taken  every  0.04  second, 
which  implies  that  40  cycles  of  a  1,000  Hertz  modulation  frequency  will  be  sampled  in 
each  interval.  This  should  give  a  good  P-P  signal  value.  The  bandwidth  of  the  filter  circuits 
that  isolate  signals  on  the  output  of  the  PMT  will  then  be 


1  1 

BW  =  -  =  -  =  12.5  Hz  (40) 

2  r  2(1/25) 


This  should  give  good  noise  rejection. 


A  scan  at  one  angular  orientation  of  the  object  transparency  will  take  33  turns, 
or  33  seconds,  to  go  from  0  to  100  cycles  per  millimeter  in  spatial  frequency.  Allowing 
a  turn-around  time,  it  is  reasonable  to  assign  40  seconds  to  a  single  radial  scan  in  spatial 
frequency. 

In  addition  to  scanning  in  spatial  frequency,  the  device  must  scan  in  angle  <p  in 
order  to  read  out  the  two-dimensional  Fourier  transform  of  the  object.  Because  of  the 
symmetry  of  the  modulus  of  the  Fourier  transform  of  a  real  object,  only  a  scan  over 
an  angle  <t>  =.  180*  need  be  taken  if  the  modulus  of  the  transform  is  of  interest.  To 

calculate  the  increment  Ax  needed  to  satisfy  Rayleigh’s  criterion,  notice  that  the  number 
of  radial  samples  just  calculated  is  given  by 


100  cycles/mm 

-  =  820  samples 

0.122  cycles/mm 


Then,  around  the  circumference  of  the  transform  plane  from  0  to  180  *  we  have 


8207r^2580  resolvable  A <t>  samples 


29 


then, 


180“ 

A0  =  - -  =  0.0698“  =  4  11  " 

2580 


The  total  number  of  data  points  will  be 


— (820)2  =  1,056,203 


and  the  total  time  required  will  be 


40  seconds  x  2580  =  103,200  seconds 


=  29  hours 


By  using  a  small  sampling  aperture,  the  number  of  sample  points  will  be  reduced  in  pro¬ 
portion  to  its  area.  Although  it  might  be  possible  to  go  to  faster  data  sampling  rates, 
this  system  has  a  signal-to-noise  ratio  problem;  therefore,  a  conservative  design  is  de¬ 
sirable.  Implementating  a  spiral  scan  through  the  transform  space  may  enable  data  to  be 
collected  more  efficiently  by  eliminating  the  redundancy  at  low  spatial  frequencies, 
which  the  radial  scan  provides.  Radial  scanning  is  assumed  in  this  design  for  experimental 
simplicity.  Current  plans  for  data  handling  are  to  plot  sine  and  cosine  transforms  on  an 
analog  plotter.  In  a  later  generation  device,  these  transforms  can  be  digitized  and  the 
Fourier  transform  phase  and  modulus  can  be  calculated  using  the  mathematical  equiv¬ 
alent  of  an  Argand  diagram.  An  Argand  diagram  is  shown  on  the  next  page. 
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DISCUSSION 


The  FTA  device  is  unique  in  its  capability  to  take  a  complete  Fourier  transform  (am¬ 
plitude  and  phase)  of  the  intensity  transmittance  of  an  input  transparency  by  purely 
optical  means.  This  capability  is  valuable  because  of  the  deleterious  effect  of  phase  noise 
on  readings  obtained  with  OPS  systems,  which  detect  the  modulus  of  the  Fourier  trans¬ 
form  of  the  amplitude  transmittance  of  the  input  transparency.  The  fact  that  the  FTA 
device  can  detect  the  entire  Fourier  transform  and  not  just  the  modulus,  may  be  signifi¬ 
cant  if  it  is  desired  to  enter  FTA  data  into  a  computer  for  further  processing,  but  it  is 
not  significant  for  pattern  recognition  work  where  the  position-invariance  property  of 
the  transform  modulus  is  a  desirable  characteristic.  As  shown  in  the  Analysis  section, 
the  ability  to  detect  a  complete  Fourier  transform  also  enables  the  subtraction  of  the 
transform  of  the  sampling  aperture  that  is  convolved  with  the  d.c.  spike  in  the  transform 
plane.  This  capability  can  be  useful  for  those  cases  where  measurements  of  low  spatial 
frequency  components  are  impeded  by  the  ringing  of  the  sampling  aperture  transform. 

Although  it  has  the  above  stated  advantages,  the  FTA  has  inherent  disadvantages, 
which  cause  problems  that  do  not  have  solutions  at  this  time. 

Briefly,  these  problems  are: 


1.  The  FTA  is  not  a  parallel  processor  in  that  the  Fourier  transform  of  the  input 
must  be  built  up  one  point  at  a  time. 

2.  Data  rates  are  slow.  This  can  be  improved  somewhat  by  engineering,  but  will 
remain  a  problem  because  of  the  mechanical  motions  that  are  required. 

3.  The  device  output  is  double  valued.  Aperturing  in  the  transform  plane  can  eli¬ 
minate  this  problem  at  a  cost  of  signal-to-noise  ratio  as  shown  in  the  Object 
Rotation  Effects  section. 

4.  The  signal-to-noise  ratio  is  poor  because  of  the  large  background  bias  that  is 
always  present. 

5.  Extremely  precise  mechanical  movements  are  required  for  reasonable  sampling 
conditions  as  shown  in  the  preceding  section.  This  requires  a  great  deal  of  stability 
in  the  apparatus. 


In  bench  work  with  a  prototype,  problems  4  and  5  were  found  to  be  very  severe. 
Although  this  is  to  be  expected  from  new  experiments,  the  degree  of  rigidity  and  pre¬ 
cision  required  to  make  repeatable  readings  would  necessitate  extreme  measures  in  a 
final  device  design. 


This  study  has  been  terminated  because  of  the  problems  listed  above.  The  insights 
into  the  details  of  Fourier  transformation  by  optical  systems  that  have  been  gained  in 
this  study  have  been  applied  in  other  projects  at  ETL.  In  particular,  a  deeper  appreciation 
of  the  problem  of  phase  noise  in  OPS  systems  has  been  achieved  that  will  lead  to  other 
less  cumbersome  methods  for  eliminating  this  noise  from  measurements  taken  with  OPS 
systems.  A  significant  reduction  or  elimination  of  phase  noise  should  lead  to  a  much 
more  valuable  pattern  classifier  that  is  capable  of  a  wider  range  of  classification  than  is 
possible  with  current  systems. 


